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ABSTRACT 


We present high resolution 11.7 and 18.3 pm mid-IR images of SN 1987A obtained 
on day 6526 since the explosion with the Thermal-Region Camera and Spectrograph 
(T-ReCS) attached to the Gemini South 8m telescope. The 11.7 pm flux has increased 
significantly since our last observations on day 6067, The images clearly show that all 
the emission arises from the equatorial ring (ER). Nearly contemporaneous spectra 
obtained on day 6184 with the MIPS at 24 gm, on day 6130 with the IRAG in 3.6- 
8 gm region, and on day 6190 with the IRS in the 12-37 pm instruments on board 
the Spitzer Space Telescope’s show that the emission consists of thermal emission 
from silicate dust that condensed out in the red giant wind of the progenitor star. 
The dust temperature is lOOii® K, an d the emitting dust mass is (2.6l 2 ;4) x 
Mq. Lines of [Ne II] 12.82 gm and [Ne III] 15.56 gm are clearly present in the Spitzer 
spectrum, as well as a weak [Si II] 34.8 pm line. We also detect two lines near 26 pm 
which we tentatively ascribe to [Fe II] 25.99 gm and [0 IV] 25.91 pm. Comparison 
of the mid-IR Gemini 11.7 gm image with X-ray images obtained by Chandra, UV- 
optical images obtained by HST, and radio synchrotron images obtained by the ATC A 
show generally good correlation of the images across all wavelengths. Because of the 
limited resolution of the mid-IR images we cannot uniquely determine the location 
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or heating mechanism of the dust giving rise to the emission. The dust could be 
collisionally heated by the X-ray emitting plasma, providing a unique diagnostic 
of plasma conditions. Alternatively, the dust could be radiatively heated in the 
dense UV-optical knots that are overrun by the advancing supernova blast wave. 
In either case the dust-to-gas mass ratio in the circumstellar medium around the 
supernova is significantly lower than that in the general interstellar medium of the 
LMC, suggesting either a low condensation efficiency in the wind of the progenitor 
star, or the efficient destruction of the dust by the SN blast wave. Overall, we 
are witnessing the interaction of the SN blast wave with its surrounding medium, 
creating an environment that is rapidly evolving at all wavelengths. Continuous 
multiwavelength observations of SN 1987A such as these provide unique snapshots of 
the very early evolution of supernova remnants. 

Subject headings: Stars: Supernovae: Individual: SN 1987A — Infrared: ISM: Dust, 
Supernova Remnants 


1. INTRODUCTION 

Since its explosion, SN 1987A has evolved from a supernova (SN) dominated by the emis- 
sion from the radioactive decay of 56 Co, 57 Co and 44 Ti in the ejecta to a supernova remnant 
whose emission is dominated by the interaction of the supernova blast wave with its surround- 
ing medium. The medium surrounding the SN is dominated by the well-known “circumstellar 
envelope” (CSE), which consists of an inner equatorial ring (ER) flanked by two outer rings (?), 
possibly part of an hour-glass structure. 

The collision between the ejecta of SN 1987A and the ER predicted to occur sometime in 
the interval 1995-2007 (??) is now underway. At UV-optical (UVO) wavelengths, “hot spots” 
have appeared inside the ER (?), and their brightness varies on time scales of a few months (?). 
New hot spots continue to appear as the whole inner rim of the ER lights up. The visible-light 
HST image obtained in 2004 Feb 20 reveals a necklace of such hot spots which nearly fill a lighted 
ring. Ongoing monitoring at X-rays with the Chandra and at radio frequencies, shows that the 
evolution of the emission from the ER follows a similar pattern at all wavelengths. 

There exist very few mid-infrared (IR) observations of supernovae in general. Therefore 
SN 1987A, the closest known supernova in 400 years, gives us an opportunity to explore the mid- 
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IR properties of supernovae and the dust in their ejecta and surrounding medium with the help of 
the newest generation of large-aperture telescopes and sensitive mid-IR instrumentation like the 
T-ReCS on the Gemini, in combination with IR data obtained from Spitzer Space Telescope(SST) 
(?). The T-ReCS observations of the mid-IR emission from SNR 1987A axe part of our continuous 
monitoring of the SN and its surrounding medium. The first detection and analysis of mid-IR 
emission at the position of the supernova has been reported in (?) (hereafter Paper I). 

The origin of the mid-IR emission could be line emission from atomic species, synchrotron 
or free-free continua, or thermal emission from dust which is probably the dominant source of 
emission. In general, there are several scenarios for the origin and the heating mechanism of 
the dust giving rise to the late time mid-IR emission in Type II supernovae (??). Thermal mid- 
IR emission could be: (1) the emission from SN-condensed dust that is collisionally heated by 
reverse shocks traveling through the SN ejecta; (2) the emission from circumstellar /interstellar 
dust heated by the interaction of the expanding SN blast wave with the ambient medium; and 
(3) the delayed emission (echo) from circumstellar dust radiatively heated by the early UVO 
supernova light curve. 

Our imaging observations provide strong constraints on these possible scenarios for the mid- 
IR emission. They show that the bulk of the mid-IR emission is not concentrated on the center of 
the explosion, but arises from a ring around the SN. The morphology of the mid-IR emission can 
therefore be used to eliminate several scenarios for its origin. Combined with Spitzer spectroscopy, 
our observations can be used to determine the dust composition and temperature distribution in 
the ring. The Chandra X-ray, and HST \JVO data provide important constraints on the physical 
conditions of the medium overtaken by the supernova blast wave. These constraints can be used 
to determine the physical association and the heating mechanism of the dust giving rise to the 
mid-IR emission. 

The paper is organized as follows: We first describe in §2 the imaging observations of 
the SN obtained by T-ReCS, and compare its IR morphology to that at radio, X-ray, and 
UVO wavelengths. Lower resolution Spitzer mid-IR imaging observations detected SNR 1987A 
(including the ring) as an unresolved point source, and are used in conjunction with spectroscopic 
data to determine the possible contribution of lines, and the composition of the dust giving rise to 
the continuum emission. In §3 we describe the procedure used for the analysis of the spectroscopic 
and imaging data, determining the dust composition, and presenting maps of dust temperature, 
IR opacities, and dust column densities. The IR image of the circumstellar medium around 
SN 1987A has a morphology similar to the Chandra X-ray and HST UVO images. The limited 
mid-IR resolution does not allow us to unambiguously determine whether the dust resides in 
the X-ray emitting gas, or in the UV-optical line emitting knots in the ER. We therefore resort 
in §4 to an analysis of possible dust heating mechanisms and to calculating the inferred dust 
masses and dust-to-gas mass ratios for several possible scenarios. In §5 we discuss the evolution 
of the supernova and its environment as manifested from the observed light curves at various 
wavelengths. The results of our paper are summarized in §6. 



- 4 - 


2. OBSERVATIONS 
2.1. Mid-Infrared Gemini Observations 

The T-ReCS mid-IR imager/spectrometer at the Gemini 8m telescope offers a combined 
telescope and instrument with diffraction limited imaging (~ 0.3” resolution) and superbly low 
thermal emissivity. On 2003 Oct 4 (day 6067), we imaged SNR 1987A with T-ReCS as part 
of the instrument’s System Verification program and we reported on the detection of 10 and 
20 pm emission from the ER, and on a 10 p m emission from the supernova’s ejecta (Paper I). 
Subsequent observations were carried out in January 6, 2005 (day 6526) in the narrow Si5 filter 
(A e ff = 11.66 p m; > 50% transmission at A = 11.09-12.22 pm), and in February 1, 2005 (day 
6552) in the Qa filter (A e // = 18.30 pm; > 50% transmission at A = 17.57-19.08 pm). Results are 
presented in Figures 1-a and 1-b. These images show several luminous “hot” spots distributed 
over the ring. The calibrated flux density integrated within an aperture of 1.3 arcsec radius 
is F 1/ (11.7 pm) — 18.4 ±1.2 mJy in the Si5 filter, and F„(18.3 pm) = 53.4 ±9 mJy in the Qa 
filter. No color correction was applied and this would most likely increase the flux density. The 
standard star used for the calibration of the 11.7 pm measurement was HD 29291, whose flux 
density was taken to be 6.78 Jy at 11.7 pm. We used a CMa with a flux density of 44.3 Jy at 
18.30 pm for the flux calibration of the 18.3 pm observation. 

The black body colour temperature corresponding to the measured fluxes at these two 
wavelengths is T = 185 K and the luminosity is Lbb = 3.74 x 10 36 erg s -1 . We use the ? 
extinction law with Tis.z = Di. 7/1-35 and Ais.zfAj = 0.083 and An^fAj — 0.098 to compute the 
black body temperature and the optical depth for each individual pixel, resulting in the maps 
shown in Fig 1-c and 1-d. Note that in order to have the iterative algorithm converge we must 
assume reasonable values as a starting point. We stress that Fig 1-c and 1-d show the color 
temperature and optical depth maps, which are slightly different from the maps related to the 
physical dust as calculated in Section 3. They are shown only for illustrating the results obtained 
from our data fitted to the simplest modeling (eg. black body). 

In Figure ?? we compare our 11.7 pm new data with the one obtained in the broad N 
band (10 pm) on October 2003 at day 6067 (Paper I). This figure shows a clear brightening 
in the South-West region of the ER, superimposed on a general brightening all over the ring. 
Given that the two images are taken with significantly different filters we investigated whether the 
difference could be due to the different spectral coverage: the spectra obtained with Spitzer which 
are discussed in next section do not show any feature that could be a source of the difference. 
Thus we conclude that we are observing a true brightness change. 

Figure ??(a,b) displays our images in both filters with the contours of the 0.3-8 keV X-ray 
image from the Chandra Observatory obtained nearly simultaneously at day 6529 (January 9-13, 
2005) (?). The correspondences between our 11.7 pm and that obtained by the Chandra is very 
good, but less so at 18.3 pm because of the lower signal-to-noise ratio in that image. A more 



-5- 


detailed discussion of the relation between the 1R and X-ray images of the ER will be presented 
in §3 below. 

Figure ??(a, b) shows the contours of the image obtained in the 12 mm band (16-26 GHz) 
at day 6003 (July 31, 2003) with the Australian Telescope Compact Array (ATCA) at the Aus- 
tralian National Telescope Facility (ATNF) (?) superimposed on our 11.7 and 18.3 pim T-ReCS 
images, respectively. The correspondences between our 11, .7 ^m image with the synchrotron 
radio emission is not as remarkable as it is with the Chandra image. Nevertheless, the 18.3 jam 
bright spot in the East side looks better correlated with the radio lobe than X-ray spots do. 
Furthermore, an image obtained in May 5, 2004 (day 6298) at the same frequencies is posted 
on the ATNF web page and is reproduced here in Figure ??(c,d). This later radio image shows 
better correlation with our 11.7 jam image than the radio image obtained on day 6003, probably 
because it was taken closer to the epoch of the mid-IR observations. This demonstrates the 
importance of evolutionary effects on the morphology of the emission at all wavelengths. 

There is good overall agreement in shape and size between our IR image and images obtained 
from the X-ray to the radio. The mean radii and approximate surface brightness distribution 
(brighter on the east side) of the ring are similar at all wavelengths, demonstrating that the dust 
is co-extensive with the gas components. The origin of that brightness asymmetry may be related 
to an asymmetric distribution of the ejecta or the CSM (??), and/or to a time-dependence effect 
caused by the tilt of the ER as argued by ?. 

The most likely source for mid-IR radiation is thermal emission from warm dust (see dis- 
cussion below). The X-ray radiation is thermal emission from very hot gas (?, and previous 
references therein), whereas the optical emission arises from the dense knots in the ER that are 
overrun by slower shocks. The radio emission is likely to be synchrotron radiation from shock- 
accelerated electrons spiralling in the remnant’s magnetic field as stated by ? and ?. ?? argue 
that, until 2000 December, hard (E > 1.2 keV) X-ray and radio emissions were produced by fast 
shocks in the CS HII region while the optical and soft (E < 1.2 keV) X-ray emissions came from 
slower shocks in the denser ER. ? note that as of 2002 Dec. 31 (day 5791) correlations between 
the X-ray and the optical/radio images are more complex than the above simple picture, which 
is expected as the blast wave is reaching the main body of the inner ring. 


2,2. Spitzer Observations 
2.2.1. Imaging Data 

Imaging of SNR 1987A was carried out with Spitzer Space Telescope’s MIPS instrument at 
24 jam (?) and IRAC instrument at 3.6-8 jam (?) (AORIDs = 5031424 and 5030912). Almost 
one year later SNR 1987A was again imaged with IRAC, serendipitously near the edge of the field 
of observations targeting other sources (AORIDs = 11191808 and 11526400). All these data were 
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obtained from the Spitzer data archive. The corresponding images are shown in Figure ??(a-f), 
which shows also a near-IR image from the Hubble Space Telescope for comparison. At 24 and 
8 pm SNR 1987A was detected as an unresolved point source amid a field of complex cirrus 
emission. At 5.8 pm, the source appears very slightly distorted, and at 4.5 and 3.6 pm the 
SN appears to be swamped by the emission of companion stars 2 and 3. The flux densities at 
24, 8, and 5.8 pm were measured using SExtractor (?) to perform aperture photometry on the 
post-BCD images. Aperture radii used were 4, 5, and 6 pixels (4.8", 6.0", and 14.7") at 5.8, 8.0, 
and 24 pm respectively. Approximate aperture corrections of 1.10, 1.07, and 1.14 were applied 
using information from the IRAC and MIPS Data Handbooks. The calibrated flux densities are 
given in Tab. ??. 


2.2.2. Spectroscopic Data 

Observations of SNR 1987A were also performed with Spitzer 1 s Infrared Spectrograph (IRS) 
in its short (wavelength) - low (resolution), short-high, and long-high modes (?). These data were 
also obtained from the Spitzer data archive (AORID = 5031168). Figure ?? shows the slit posi- 
tions for the different sets of observations. For each spectral mode we extracted SN 1987A spectra 
from the 2-D coadded post-BCD images using SPICE (http://sse.spitzer.caltech.edu/postbcd/spice.html). 
For the short-low data, the source is placed at 4 positions along the slit to generate spectra at 
two positions for two different spectral orders. For each slit position, a 2-D background image 
was generated from the median value of the images for the other three slit positions. These 
backgrounds were subtracted prior to extracting spectra with SPICE. For the high resolution 
observations, only two slit positions are observed within a much narrower slit. So for these 
observations, SPICE was used to extract the SN spectrum from the columns occupied by the 
unresolved source, and a background spectrum from columns toward the opposite end of the slit. 

The background spectra were then subtracted from the source spectra. For all modes the spectra 
from the two slit positions (per spectral order) were combined using a weighted (l/c 2 ) average, 
and (generally noisy) data where spectral orders overlapped were discarded. For the short-hi 
data, an empirically determined scaling factor of 1.38 was applied before averaging to bring re- 
sults from both slit positions into agreement. An additional scale factors of 1.46 was applied to 
the short-high spectrum to normalize it to the short-low data, and a subsequent factor of 1.25 
was applied to long-high data to normalize to the short-high data. These latter scaling factors 
are to be expected if the subtracted background was partially contaminatead by the source, or 
in the case that the smaller high resolution slits were slightly misaligned. Finally for the sake 
of comparison with the broad band measurements and dust models, the high resolution spectra 
were median-binned using intervals of 12 and 11 wavelength samples for the short and long 
wavelengths respectively. 

Figure ?? shows the overall calibrated spectrum and Figure ?? displays individual profiles 
and identifications of the main emission lines detected with IRS from the Spitzer Space Telescope: 
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it shows that the T-ReCS observations are dominated by the dust continuum emission, and not 
the lines. The [Ne II] 12.81 ja m and [Ne III] 15.56 /zm lines are however clearly seen. A weak 
[Si II] 34.81 /zm line remains after background subtraction, while [S III] 33.48 /zm disappears 
entirely. Two strong lines are seen near 26 /zm: the redder line is [Fe II] at 25.99 /zm and the 
bluer line is [0 IV] at 25.89 /zm. Both lines could be arising from [Fe II] (or both [0 IV]) if from 
fast moving ejecta on near and far sides of the explosion, but the lack of splitting of the other 
lines makes this seem unlikely. A weak third component seems present here as well. If fitted as 
an unresolved line this may be [F IV] at 25.83 /zm. However, it may also be fitted as the wing 
of a very broad line (FWHM ~ 2800 km s _1 ) underlying the narrower [O IV] and [Fe II] lines. 
Such high velocity would indicate an association with the SN ejecta rather than the ER for this 
line. Line fluxes, centers, and widths have been calculated by gaussian fits using SMART (?) 
and results are given in Tab. ??. 


2.3. Near-IR CTIO Observations 

Near-IR J(1.25 /zm), H(1.65 /zm), and K(2.2 /zm) imaging observations of SNR 1987A were 
obtained on 2005, January 3-5 with ISPI attached to the Blanco 4-m telescope at the Cerro 
Tololo Interamericain Observatory. Results are displayed in Figure ??. We clearly see the ER at 
J and K, while most of the emission detected in the H band arises from the supernova itself. These 
images have been deconvolved by using the “Multiscale Maximum Entropy Method” algorithm 
developed by ?. The flux densities measured in the CTIO data are ( Fj = 0.33 ± 0.25 mJy; 
Fjj = 0.11 ± 0.05 mJy; Fk — 0.51 ± 0.2 mJy). These fluxes refer to the total of the ER + ejecta 
in each band (although each band seems dominated by one or the other). Images obtained with 
HST on November 11, 2005 through the filters F110W, F160W, and F205W, are also shown in 
this figure, and show that compared to the ring, the ejecta are relatively brighter in H band 
(F160W) than in other bands. This clearly demonstrates the validity and the power of the 
deconvolution algorithm applied to our CTIO data. Note, however, that the HST bands are not 
exactly J, H, K: one significant difference is that the Paschen a line is missed between the ground 
based H and K bands, but would be seen in HST’s F205W band. 

The Spitzer mid-IR spectrum indicates that dust in SN 1987A and the ER is too cool (see 
Section 3) to emit significantly in these near-IR bands. We believe that line emission dominates 
these broadband near-IR observations, and that differences in the composition, temperature, and 
density of the ejecta vs. CSM lead to the different relative brightnesses of these structures. The 
young supernova remnants Cas A and Kepler may be the closest analogs we have for interpreting 
SN 1987A’s near-IR emission. Near-IR spectra and imaging of these SNRs by ? and ? reveal 
similar variation in near-IR colors. Circumstellar material (e.g. Cas A’s quasi-stationary flocculi) 
is bright in J due to He 1 1.083 /zm, and emits strong [Fe II] lines in both J and H bands. Weaker 
hydrogen lines can be detected in all bands. The ejecta in several of Cas A’s fast-moving knots 
emit strong [S II] lines at 1.03 jam in J, and weaker lines in [Fe II] and Si in various ionization 
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states. Thus, the near-IR emission of the circumstellar ER may be dominated by H and He I lines 
in J and H with synchrotron emission possibly augmenting the K band (??). The synchrotron 
processes deserve consideration in SN 1987A particularly in view of the fact shown in Figure ?? 
that the radio morphology of SN 1987A reproduces the ring shape now. Also, the H2 molecule 
has many emission lines in the K-band and longer wavelengths. Spectra are clearly needed for 
testing the plausibility of the contribution from this molecule or otherwise. The first overtone 
band head of molecular CO lies at 2.29pm and may be also contributing to K-band emission 
in the ring. The relatively strong H band emission of SN 1987A’s ejecta could reflect [Fe II] 
emission, because a large fraction of the ejecta should be Fe from the core of the explosion, 
whereas Cas A’s FMKs generally consist of lighter metals originating in the outer layers of the 
progenitor. 

3. DATA ANALYSIS 
3.1. The Dust Properties 

Figure 5 clearly demonstrates that the mid-IR emission is dominated by thermal emission 
from dust. The specific luminosity of a single dust particle of radius a at temperature T d , at 
wavelength A is given by: 


4(A) = 47ca 2 7rB„(\, T d )Q{ A) 

= 4m d /c(A, T d ) 


( 1 ) 


where B„( A, T d ) is the Planck function, Q( A) the dust emissivity at wavelength A, and «( A, a) = 
3Q(A) / 4pa is the dust mass absorption coefficient, where p is the mass density of the dust particle. 
In the Rayleigh limit, when a < A, k is independent of particle radius. Figure ?? illustrates the 
k( A) curves for amorphous carbon (?), graphite and silicate grains (?) over the 5 to 30 pm 
wavelength region. The figure shows the distinct optical properties of these dust particles over 
the mid-IR wavelength regime which can greatly facilitate the identification of the emitting 
material with even limited broad band filters. 

For an optically thin point source, the flux density, F u ( A), at wavelength A is given by: 


F v { A) =4 M d 


k(A) ?rR v (A, T d ) 
47tD 2 


( 2 ) 


where M d is the dust mass, and D is the distance to the supernova, taken to be D = 51.4 kpc 

(?)• 

For an extended optically thin source with an angular size 0, the surface brightness, J„(A), 
is given by: 


4(A) = Q r d (X)B u (X, T d ) 


( 3 ) 
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where r d ( A), the dust optical depth, is given by: 


T d (X) = 


M d k { X) 
SID 2 


(4) 


3.2, Spectral Analysis 

We fitted the integrated T-ReCS flux densities with a population of dust particles consisting 
of a single population of bare graphite, silicate or amorphous carbon grains, using equation (2). 
Optical properties of silicate and graphite grains were taken from ?, and those for amorphous 
carbon (BE) were taken from ?. The results are given in Figure ??, which shows that the T-ReCS 
observations alone cannot discriminate between the different dust compositions. Comparison 
with the Spitzer observations clearly demonstrates that the IRS observations can be well fit 
with a silicate dust composition (Figure ??), ruling out graphite or carbon dust as major dust 
constituents in the CSM. 

Figure ?? shows the fit of the silicate dust spectrum to the Spitzer IRS spectrum. The 
silicate temperature is T d — 1801^ K, and the dust mass is M d = (1. lioJ) x 10 -6 M©. The global 
parameters resulting from our model for silicates are given in Table ??. Figure ?? shows the 
residuals of the fit, obtained by subtracting the silicate model fit from the data. The residuals 
are small, and their spectrum is too sharply peaked at short wavelengths to be fitted with any 
blackbody. The residuals are also broader than typical atomic line widths, suggesting that they 
may be due to solid state features in the dust, reflecting differences in the crystalline structure of 
the silicates in the CSM from the average interstellar silicate dust used in the model. This figure 
shows also that the residuals cannot be fit by emission from amorphous carbon dust (shown by 
the green line in this figure) . 


3.3, Image Analysis 

Our spectral analysis showed that the mid-IR spectrum of the ER is dominated by silicate 
emission. We therefore used the 11.7 and 18.3 pm. images of the remnant to construct tempera- 
ture, and dust opacity maps of the circumstellar ring using Eq. 3 to derive the dust temperature, 
and Eq. 4 to derive the dust optical depth. In calculating these quantities, we applied a back- 
ground threshold of 0.03 mJy/pix at 11.7 pm and of 0.08 mJy/pix at 18.3 pm. The mass of the 
ring was calculated from the optical depth, over a surface area of 269 pixels, corresponding to 
the number of pixels that had a flux exceeding the respective thresholds at each wavelength. The 
average dust temperature in these maps is T d = 1661}® K. The average 11.7 pm optical depth 
per pixel is r d = (5.5 j^f) x 10 ' 6 , and the total dust luminosity is L d = (2.3lto!4) x 10 36 erg s _1 , 
giving a dust mass M d = (2.61^4) x 40~ 6 M©, in good agreement with the total mass obtained 
from the spectral analysis of the ER. Figure ?? shows the maps of the silicate dust temperature 
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(a) and optical depth (b) in the ER. 


4. THE ORIGIN OF THE MID-INFRARED EMISSION 

The data presented in this paper show unambiguously that the emission is thermal emission 
from dust. At issue are the location and heating mechanism of the dust. The forward expanding 
non-radiative blast wave is currently interacting with the circumstellar material and the knots 
in the ER. The interaction of the blast wave with the knots transmits lower velocity radiative 
shocks into these dense regions, producing soft X-rays and the “hot spots” seen in the HST 
images. The interaction of the blast wave with the dense knots also generates reflected shocks 
that propagate back into the medium that was previously shocked by the expanding SN blast 
wave. The complex morphology and density structure of the ER gives rise to a multitude of 
shocks characterized by different velocity, temperatures, and post shock densities. 

The mid-IR images cannot determine the location of the radiating dust, whether it resides 
in the X-ray emitting gas or in the denser UVO emitting knots. Therefore, we can not, a priori, 
assume a particular dust heating mechanism: collisional heating in the shocked gas, or radiative 
heating in the radiative shocks. 


4.1. Dust Heating Mechanism 

The relative importance of the two dust heating mechanisms is given by the ratio (?): 

rn Hrad U. e fljj A(T' e ) Pi ibs 

= H^i~ n~d n e A d (T e ) 1 j 

where H ra d and H co u are, respectively, the radiative and collisional heating rates of the dust, 
n e , T e , and A(T e ) are, respectively, the electron density, temperature, and the atomic cooling 
function (erg cm 3 s -1 ) of the gas, P a bs is the fraction of the cooling radiation that is absorbed 
by the dust, and A d(T e ) is the cooling function of the gas via electronic collisions with the dust 
and given by: 

A d {T e ) = 2 v e 7ra 2 kT e {h) (6) 

where rid the number density of dust grains, a their average radius, v e — (SkTe/irme) 1 / 2 is the 
mean thermal speed of the electrons, and (h) < 1 is the collisional heating efficiency of the dust 
which measures the fractional energy of the electrons that is deposited in the dust (?). 

In the following we will examine possible locations of the dust giving rise to the IR emission. 
For each possible site, the X-ray emitting gas or the optical knot, we determine the dominant 
cooling mechanism, the temperature of the dust, and the inferred dust mass. A site is viable 
if it can maintain a dust temperature between ~ 150 and 200 K, the range of values reflecting 
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the range of dust temperatures derived from the T-ReCS and Spitzer observations, and if the 
derived dust mass does not violate any reasonable abundances constraints. 


4.2. Dust in the X-ray Emitting Gas 


The morphological similarities between the 11.7 pm emission, dust temperature, and optical 
depth maps on one hand and the X-ray maps of the supernova on the other hand suggests that 
the dust giving rise to the IR emission may be well mixed with the X-ray emitting gas. 


For an optically thin plasma P a bs ~ Td = rid^a 2 {Q)E , where (Q) is the radiative absorption 
efficiency of the dust averaged over grain sizes and the radiation spectrum, and i is a typical 
dimension of the emitting region. Inserting this value in eq. (5) we get: 


n 


n e A(r e )i (q) 

2 v e kT e (h) 

n e A(T e ) i 


7r m e 
~32 


(kT e fP 



( 7 ) 


For average conditions in the X-ray plasma (see below), characterized by plasma tempera- 
tures and densities of T e .« 10 7 K, and n e « 300 — 10 3 cm -3 , we get A (T e ) «4x 10~ 23 erg cm 3 s _1 , 
and (h) > 0.1 for dust particles with radii larger than 0.05 pm. Adopting a value of (Q) ~ 1, 
gives an upper limit of 

TZ<2.1 x 10 _2 ° l (8) 

The size of the X-ray emitting plasma is less than the radius of the ER which is i < 0.7 lyr 
= 7 x 10 17 cm, giving 7 Z < 0.14, that is, the dust heating in the X-ray gas is dominated by 
electronic collisions. 


4-2.1. The Dust Temperature 

Figure ?? depicts contours of the temperature of collisionally heated silicate grains of radius 
a — 0.10 pm as a function of plasma density and temperature. To calculate the energy deposited 
in the dust we used the electron ranges of ? for energies between 20 eV and 10 keV, and those of 
? for higher incident electron energies. We also assumed that all incident electrons penetrate the 
dust (no reflection). Above gas temperatures of about 3 x 10 6 K, the temperature of collisionally 
heated dust is independent of grain radius, and very well represented by a single dust temperature, 
since both the radiative cooling and the collisional heating rates are proportional to the mass 
of the radiating dust particle. Furthermore, at these plasma temperatures the dust temperature 
is essentially independent of gas temperature as well, and therefore an excellent diagnostic of 
plasma densities. The figure shows that for plasma temperatures above ■'v 3 x 10 6 K, dust 
temperatures between 150 and 200 K require electron densities of about 300 to 1400 cm -3 . 
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X-ray observations show the presence of two X-ray emission components (?) : one associated 
with the “slow shock” with an electron temperature of ~ 0.23 keV and a density of n e ~ 
6000 cm -3 , and the second, associated with the “fast shock” with electron temperatures and 
densities of kT ~ 2.2 keV and n e ~ 280 cm -3 , respectively. The lastest Chandra data indicate 
(?) that the “fast shock” and the “slow shock” of the model are becoming less distinguishable, 
as the overall shock front is now entering the main body of the inner ring, that is, the electron 
temperature of the soft component is increasing and that of the hard component is decreasing. 
Albeit rather speculative, the overall temperature might thus be “merging” onto an “average” 
temperature of kT ~ 1.5 keV ~ 1.8 x 10 7 K, ?, and intermediate electron densities. The range 
of densities expected for this “average” shock is in very good agreement with that implied from 
the IR observations. 


4-2.2. The Infrared-to-X-ray Flux Ratio (IRX) 

An important diagnostic of a dusty plasma is the infrared-to-X-ray (IRX) flux ratio (?). 
For a given dust-to-gas mass ratio and grain size distribution, the IRX ratio is defined as the 
infrared cooling to X-ray cooling in the 0.2-4.0 keV band, and is given by: 


IRX(T e ) 


Tie A d(T e ) 
n.# A x (T e ) 

irninZd Ad(T) 

(m d ) A x (T e ) 


( 9 ) 


were Z d = nd(md) / pniinriH is the dust-to-gas mass ratio, p, is the mean atomic weight of the gas, 
and ( md ) the mass of a dust particle, averaged over the grain size distribution. For a given value 
of Z d , the IRX ratio is only a function of plasma temperature, and ranges from a value of about 
10 for T e = 10 6 K to a value of ~ 400 for T e = 10 s K. For young supernova remnants, ? found 
that the IRX ratio is significantly larger than unity for 7 of the 9 remnants considered in their 
paper, the other two having only an upper limit on their IR emission. 


The observed IRX ratio can be obtained from the X-ray and IR fluxes from the SN. From 
the January 2005 Chandra data, we estimate the X-ray flux in the 0. 2-4.0 keV X-ray band to be 
7.24 x 10~ 12 erg cm -2 s _1 after correcting for interstellar absorption by an H-column density of 
N h = 2.35 x 10 21 cm- 2 . 

Note that these values are based on the two-temperature model as used in ?. Fractional 
contributions from the soft (kT e — 0.3 keV) component in the total flux is ~70% for the unabsorbed 
flux (and thus for Lx)- Thus, a contribution from each component in the 0.2-4 keV X-ray flux 
seems to be significant rather than being dominated by one of them. 

The total IR flux is Fir — 7.7 x 10 -12 erg cm“ 2 s" 1 which leads to an IRX ratio ~ 1. This 
value is lower than the values reported in Paper I, in which IRX = 6 for the decelerated slow 
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shock component, and IRX = 3 for the blast wave shock in the two-temperature model. It is 
much lower than the theoretical value of ~ 10 2 , expected for a T e « 1.8 x 10 7 K plasma, which is 
observed in the young remnants Tycho and Cas A. Other (mostly older) SNEs with measurable 
IR emission show somewhat lower IRX ratios, but not as low as SNR 1987A. 

Several effects could be the cause for this very low value of the IRX ratio in SNR 1987A. 
First, in remnants the IRX ratio was calculated by ? for an interstellar dust-to-gas mass ratio 
of 0.0075. The SNR 1987A blast wave is expanding into the circumstellar shell of its progenitor 
star, which, a priori, is not expected to have an interstellar dust-to-gas mass ratio. Moreover, 
when estimating the depletion of elements onto dust in the ER, we should compare the expected 
dust abundances in the ER with the maximum available abundances for the LMC. General 
LMC abundances exist for B stars (?) and ISM (?) and, although controversy still remains 
(?), the LMC metallicity is usually assumed to be 0.5 — 0.7 solar. Assuming that the fraction 
of metals locked up in LMC dust is the same as in the local ISM, and that the ER has the 
same metallicity as the LMC, then the IRX ratio in the ER should be about 0.5 — 0.7 times 
that expected from Supernova Remnants in the Milky Way, still significantly larger than implied 
from the observations. The extremely low value of the IRX ratio may therefore be due to a 
deficiency in the abundance of the dust, compared to interstellar values, which may reflect the 
low condensation efficiency of the dust in the circumstellar envelope. 

Second, the dust deficiency could be the result of grain destruction by thermal sputtering 
in the hot gas. The sputtering lifetime, T sput , in a plasma with temperatures above 10 6 K is 
about (?): 

Tsput « 3 x 10 5 yr (10) 

n(c m~ 6 ) 

where n is the density of nucleons in the gas. The X-ray emitting gas is highly ionized, and we will 
assume that its density is that required to heat the dust to the observed range of temperatures, 
that is, n « 300 to 1400 cm -3 . Grain destruction is important when the sputtering lifetime is 
about equal to the age of the shocked gas, which we take to be ~1 yr. The low IRX ratio can 
therefore be attributed to the effects of grain destruction if the dust particles had initial radii 
between 10 and 50 A. So attributing the small IRX ratio to the effect of grain destruction in the 
hot plasma requires that only small grains had formed in the presupernova phase of the evolution 
of the progenitor star. 

The low IRX ratio shows that IR emission from collisionally heated dust is not the dominant 
coolant of the shocked gas. Its lower than expected value suggests a dust-to-gas mass ratio in 
the ER that is only a few percent of its interstellar value. The puzzle of the low dust abundance 
is greater if, in fact, the IR emission arises from dust that is not embedded in the X-ray emitting 
gas, but from dust that resides in the UV/optical knots instead. 
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4.3. The Dust Heating mechanism 

The possibility that the initial size of the dust grains swept up by the shock is small suggests 
that the temperature of the grains may not be at the equilibrium value but may fluctuate due 
to the stochastic nature of the heating and cooling. Very small grains will be stochastically 
heated if the energy deposited in the grain in a single collision is large compared to its internal 
energy content, and if its cooling time via IR emission is shorter than the time between subsequent 
collisions. Assuming that the electron and ion temperature is instantaneously equilibrated behind 
the shock, the mean thermal energy of the electrons will be 2.6 keV at a postshock temperature 
of 2 x 10 7 K. Using the electron ranges given by ? we get that the energy, A E, deposited in a 
dust grain of radius a(pm) and density p ( g cm -3 ) by electrons with energies E > 370 eV is given 
by: 

A-E(erg) = 5.5 x 1CT 7 (11) 

which for p = 3 g cm -3 , a = 0.0050 pm, and E = 2.6 keV gives an energy deposition A E = 
1.7 x 10 _1 ° erg. This value is small compared to the internal energy of the grain, which is about 
6 x 10 -10 erg at the equilibrium grain temperature of 170 K. The average thermal speed of the 
electrons is about 3 x 10 9 cm s -1 , giving an average time between electronic collisions of about 
1 sec for an electron density of 300 cm -3 and a grain radius of 0.0050 pm. This is shorter than 
the grain’s cooling time at 170 K, which is about 10 s (?)., suggesting that the dust can maintain 
an equilibrium temperature at that value. 

Collisions with protons can deposit a significantly larger amount of energy in the dust. 
Typical proton energies in the shocked gas are ~ 5 keV, and the stopping power for protons 
with that energy is about 240 MeV cm 2 g -1 (NIST tabulated values: http://physics.nist.gov/cgi- 
bin/Star). The energy deposited in a 0.0050 pm radius dust particle is then about 6 x 10~ 10 erg, 
which is about equal to the internal energy of the dust grain at 170 K. So collisions with protons 
and heavier nuclei can be neglected if they are in thermal equilibrium with the electrons, because 
of their lower collision rate. However, if the electron and ion temperatures are not equilibrated 
behind the shock, then the dust heating rate will be dominated by collisions with the protons. 
The time between successive collisions will be about 40 s for a proton density of 300 cm -3 . The 
grains should therefore cool to a temperature of about 100 K before another collision will take 
place. However, the data does not support such a broad range of dust temperatures, suggesting 
that the gas density could be somewhat higher (by a factor of ~ 3), in which case the dust will 
cool only to a temperature of ~ 140 K, which may be more consistent with the observations. 

All these scenarios support the idea that the dust temperature in the X-ray emitting gas 
does not fluctuate wildly about the equilibrium dust temperature, which can still provides strong 
constraints on the density and the equilibration of the electron and ion temperatures in the 
postshock gas. 
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4.4. Dust in the Dense Knots of the Equatorial Ring 

The UVO light emitting knots discovered with the HST resemble a string of beads uniformly 
distributed along the ER. Figure ?? and Figure ?? depict the map of the dust optical depth 
overlayed with contour levels of the HST emission obtained on Dec 15, 2004 (day 6502), the 
closest to our 11.7 yum observations. The data look very similar, but the IR emission seems 
to emanate from a somewhat wider region than the optical emission, an effect that cannot be 
entirely accounted for by the lower resolution of the IR data. 

Nevertheless the good correlation between the IR emission maps and the HST image, suggest 
that a significant fraction, if not most, of the mid-IR emission may be emanating from the knots. 
The physical conditions of a particular knot (Spot 1 on the ER) have been modeled in detail by ?, 
from the analysis of the U V / optical line emission detected by the HST Space telescope Imaging 
Spectrograph (STIS). They found that the UV fluxes could be fit with a model consisting of 
two shocks with velocities of v s = 135 and 250 km s -1 expanding into preshock densities of no 
— 3.3 xlO 4 and 10 4 cm -3 , respectively. The postshock temperatures behind the slow and fast 
shocks are T s = 4.5 x 10 5 , and 1.5xl0 6 K, with gas cooling rates of A (T) ~ 10 -22 erg cm 3 s -1 . 
The cooling time of the shocked gas is given by t coo i = kT s /nA(T s ), and is equal to 0.1 and 
2 yr for the slow and fast shock, respectively (?). The thickness of the shock front is therefore 
~4x 10 13 and 8xl0 14 cm for the slow and fast shocks, respectively, both significantly smaller 
than a typical radius of the knot, which is about 2 x 10 16 cm. So most of the dust in the knot 
resides in the unshocked gas and is heated by the radiation emitted from the cooling shocked 
gas. 


The radiative energy density seen by the dust is approximately given by: 


TT _ n 2 K{T s )i C ooi 

V rad ^ 

c 

« n k T s — 
c 

« fcnokTs- 
c 

where iaoi = tcooi and / c is the compression factor of the gas in the postshock 
T s — 10 6 K, no = 10 4 cm -3 , and a shock velocity of 200 km s -1 we get that 


( 12 ) 


region. For 


Urad « 2 x 10 9 x f c erg cm 3 


(13) 


? find that compression factors can be as large as « 550, giving radiation densities of ~ 
10~ 6 erg cm -3 throughout the knot. The energy density of the local interstellar radiation field is 
about 3 x 10~ 12 erg cm -3 . Silicate dust particles immersed in this field achieve equilibrium dust 
temperatures of about 15 K (?). The energy density in the knot is therefore higher by a factor 
of ~ 3 x 10 5 than that of the local interstellar radiation field, and the average dust temperature 
should therefore be higher by a factor of rsj 8.3 for a A 2 dust emissivity law (Td oc U^)'. This 
gives a typical dust temperature of ~ 125 K, in reasonable agreement with the observed average. 
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The total mass of radiating dust was found to be ~ 10~ 6 M©. The typical mass of gas in 
a knot of radius r = 2 x 10 16 cm, and density no = 10 4 cm -3 is ~ 10 -4 M©. If N « 20 is. the 
number of knots in the ER, then the dust-to-gas mass ratio is ~ 10 ~ 6 / (N x 10~ 4 ) « 5 x 10~ 4 , or 
approximately a factor of 10 less than the average dust-to-gas mass ratio in the local interstellar 
medium. 

The low abundance of dust in the knots could be explained if the dust is efficiently destroyed 
in the shocked gas and if the transmitted shocks have already traversed most of the volume of 
the knots. Calculations presented by ? show that about 49% of the silicate grains swept up 
by 200 km s _1 shocks expanding into a medium with a preshocked density of 0.25 cm -3 are 
destroyed. This fraction could be significantly higher for the densities encountered by the shocks 
traversing the knots. This scenario predicts that the IR emission from the knots was higher in 
the past, contrary to observed IR light curves (see below). Therefore, if the IR emission emanates 
from the knots, the low dust abundance must reflect the initial dust abundance in these objects. 

5. THE LIGHT CURVES 

The light curves at 10 and 20 /am are shown in Figure ??: the absolute flux calibration have 
been made using ? 0-magnitude fluxes [Fo(10.0 fx m) — 35.24 mJy, Fo(11.7 fx m) = 28.57 mJy, 
and Fo(18.30 fxm) = 10.25 mJy]. It can be seen from Figure ?? that the flux arising from the 
ejecta at 10 gm declines exponentially from day 2200 through day 4200 (the ISOCAM and OSCIR 
observations) until day 6000, at a rate of ~ 0.32 mag y~ l . This would imply that the observations 
on day 4200 are not dominated by the ring emission. It is likely that the ring emission started 
around day 4000, at roughly the epoch when the first optical spot was discovered (?), and in good 
agreement with Fig. 5 of ? which presents ATCA and Chandra/ROSAT data. The significant 
increase of the fluxes reported in the present paper compared to the previous fluxes (Fig. ??) 
is consistent with the soft X-ray flux increase observed in the last set of data (?). This clearly 
shows that “something” has happened at around day ~ 6000. 

This is also manifest in the radio light curve at 843Mhz from MOST. The earlier observations 
by (?)from 3000 to 6000 days show a steadily increasing flux together with a steadily increasing 
rate of change. Near day 6000 however (?) the rate of increase undergoes a more abrupt change 
upwards, and is therefore clearly associated with the changes seen at other wavelengths. Figure ?? 
also shows the Spitzer data point resulting from integrating the IRS data within the appropriate 
bands. Including these data points shows that the trend of the increase of the flux at 11.7 yum 
since after day ~6000 is better fit by a “linear” function, while the increase in magnitude vs. time 
is not linear. Although we do not make strong conclusions, this could indicate that the shock is 
travelling through structures with cross sections (as seen from the center of the explosion) that 
do not increase with radial distance, like short cylindrical clouds. Sputtering time scales for all 
but the smallest grains are too long (10-100 years) to affect the IR emission, and the plasma 
cooling time due to the IR emission is also too long (~ 100 years) to affect the time variation in 
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the IR emission. 

The 13 cm radio emission light curve (http://www.atnf.csiro.au/research/SN1987A/) is also 
shown in Figure ??. No data are currently available after day 6244 so it is premature to discuss 
the associated 13 cm evolution. It would be surprising if the rate of increase does not change 
upwards. While the sudden increase in the X-ray light curve at day around 3700 has been 
interpreted as the encounter of the shock front with the first protrusions of the ER .(?), the 
“jump” after day around 6000 could be the sign of the shock reaching the main body of the 
ER (?). These authors also show that the light curve of the hard X-rays (3-10 keV) is much 
flatter than one corresponding to the soft X-rays, and similar to the radio light curves, and they 
argue that it is likely that the hard X-ray emission comes from the fast reverse shock rather than 
the decelerating forward shock, just like the radio emission (?), We note that the reverse shock 
origin for the radio emission is one possible interpretation for explaining the inconsistency of the 
radio images with the IR and the soft X-rays images. It is not our intention to discuss evolution 
of radio fluxes and their relationship to other changes reported in this paper, but we note that 
the log-log plot in Figure 16 hides a significant change in the rate of increase of the 13 cm at 
5000 days. The epoch of increase in the rate appears to be frequency dependent. 


5.1. The Ejecta 

An asymmetry in the profiles of optical emission lines that appeared at day 530 showed that 
dust had condensed in the metal-rich ejecta of the supernova (?). Although it was discovered 
via spectroscopy, the presence of the dust could be easily inferred from the spectral energy 
distribution: as the dust thermalized the energy output, after day 1000 SN 1987A radiated 
mainly in the mid-infrared (?) . Although the presence of the dust emission from the condensates 
in the ejecta was reported in Paper I, we do not detect them in the present observations. The 
possibility that this was due to the fact that the present observations were achieved with the 
narrow Si5 filter instead of the broad N filter was investigated, for the occurence of the [Ne II] 
12.8 jum line. Indeed, this line corresponds to a fine structure transition in the ground state of 
Ne II so the temperature to excite the upper level does not need to be high if there is Ne II. 
However, Paper I reports a dust temperature of 90 K < T < 100 K in the ejecta, and there is 
little evidence for X-ray emission from the ejecta. Furthermore, the X-ray emitting gas would 
be ionized to a much greater extent than to produce Nell. It is thus most likely that Ne II is 
coming from warmer regions near the X-rays which are also responsible for the other HST and 
Spitzer lines. Why then is the ejecta not detected in our last observations? 

In Paper I, the N band background sigma was 0.033 mJy/pix. We used a 12 pixels area to 
integrate the central source, and then the 3-<r detection limit was 0.34 mjy which is about the 
0.32 mjy measurement reported for the central source in Paper I. In the present data, although 
the source is observed with a better signal-to-noise ratio, the background at 11.7 /xm is affected 
by a higher noise, with a 3-cr background value of 0.47 mjy. It transpires that we could not have 
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detected a source at the flux level expected from the radioactive decay based light curve (eg. 
~ 0.3 mJy). It appears then that the only detection of the ejecta at this late stage is achieved 
in the H band (see Figure ??). The origin of this emission is most likely due to line emission, 
although we note that the limiting flux for any continuum emitter at the center of SN 1987A, 
in the wavelength range 2900-9650 A reported recently by ? (4.3 mJy in the I filter) is much 
above our detection at H (0.11 mJy). Therefore, our near IR results are compatible with HST 
observations and do not exclude contribution from a continuum. 

Supernovae are known contributors to interstellar dust. The presence of isotopic anomalies in 
meterorites (?) and observations of Cas A (??) and SNR 1987A (Paper I and references therein) 
provide direct evidence for the formation of dust in SNe. However, the relative importance of SNe 
compared to quiescent outflows from AGB stars in the production of interstellar dust grains is 
still unclear (???). So far, the dust masses in the the ejecta have been determined observationally 
for two SNe, SN1987A (?) and SN1999em (?). These masses were ~ 3.10 -4 M© and ■<■*> 10 -4 M© 
. In both cases the authors note that these values could be much higher if the dust exists in 
opaque clumps which may be the case. At face value however they are much less than the 0.1 - 
1.0 M© required to make SNe the dominant source of interstellar dust particles. In future one 
must find a means of quantifying the effects of clumping. Observations of the handful of other 
supernovae in which an IR excess is interpreted as dust forming in the ejecta (SN 1979C and 
SN 1985L, and probably SN 1980K) do not allow an estimate of the mass of dust. As for the 
young Galactic SNR which have been observed by IRAS and ISO (Cas A, Kepler and Tycho) 
the dust mass deduced is only 10~ 7 — 10 -3 M© (??), also many orders of magnitude lower than 
the solar mass quantities predicted. 

Paper I discusses the role of supernovae in dust production. We showed there that the 
dust which condensed in the ejecta of SN 1987A has survived 16 years since outburst, and 
was still radiating the energy released by the radioactive decay of 44 Ti at the expected level. 
Unfortunately, we could not accurately estimate the mass of the dust, and the observations 
reported in this paper do not allow it either, in the absence of detection of the ejecta. Thus, if 
we consider SN 1987A as an archetype, our data can neither support nor rule out the hypothesis 
that supernovae are significant sources for dust production. 


6. CONCLUSIONS 

We have presented mid-IR images of SNR 1987A obtained with T-ReCS on the Gemini South 
telescope on day 6526 at 11.7 and 18.3 /u m and with IRAC (5.8 and 8 yum) on day 6130 and 
MIPS (24 yum) on day 6187 onboard Spitzer, together with 3 - 37 yum spectroscopic observations of 
the remnant obtained with IRS on day 6190 at the same observatory. The imaging observations 
(Figure ??) show that the mid-IR emission arises from the dust in the equatorial ring (ER) heated 
up by the interaction of the SN blast wave with its eircumstellar medium. Several theoretical 
models predicted the presence of dust in the CSE of SN 1987A which was produced in the winds 
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of the supergiant phase. The location of the TR emission rules out the possibility that the dust 
condensed out in the SN ejecta, strongly suggesting a circumstellar origin instead. 

The 3-37 pm spectrum (Figures ??) shows that the emission arises from a population 
of astronomical silicate particles. Temperature maps show that the dust temperature is fairly 
uniform in the Elt and about 1661 jj K, with total dust masses of ~ 2.61^4 x 10~ 6 M©. 

A comparison with Chandra and HST observations show an equally good correlation between 
the 11.7 pm IR and the X-ray and UV-optical images of the supernova. Because of the limited 
angular resolution in the IR we cannot determine the location or heating mechanism of the 
radiating dust. 

The dust could be residing in the hot ~ 10 7 K gas and collisionally heated by the X-ray 
emitting plasma. The dust temperature is then an excellent diagnostic of the electron density, 
giving a value of ~ 300-1400 cm -3 , similar to the value suggested by the Chandra observations. 

Comparison of the IR and X-ray fluxes suggests that the dust is depleted by a factor of 
~ 30 in the X-ray emitting gas, compared to its value in the local interstellar medium of the 
LMC. This low value could be due to its destruction by thermal sputtering in the shocked gas, 
requiring the initial grain radii to be below r\j 50 A. 

Alternatively, the dust could be residing in the UV-optical emitting knots in the ER and 
radiatively heated by the cooling gas that was excited by the shocks propagating through these 
knots. A simple calculation for a particular knot shows that the radiative energy density can 
heat the dust to typical temperatures of about 125 K, similar to those inferred from the IR 
observations. A comparison with the mass of the knots shows that the dust-to-gas mass ratio in 
the knots is lower by a factor of ~ 10 compared to its value in the ISM of the LMC. This low 
abundance reflects the low condensation efficiency of the dust in the outflow of the progenitor 
star. 

We stress that in order to assess the role of SNe in the production of dust in the Universe, 
it is clearly important to measure the presence of dust that survives into the formation of the 
remnant, and for this, mid-IR and sub-mm observations are critical. 
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Table 3. Global parameters from Silicate Dust Modeling 1 


Instrument 

Day 1 2 

Tdust (K) 

Ljr (erg s 1) 

Mdust (Mq) 

T-ReCS 

6070 

180 i 2 o 

9 ± 3 x 10 35 

1 - 4 x 10- 6 

Spitzer 

6190 

18011 ° 

1.6 ±0.3 x 10 36 

0.7- 1.7 x 10~ 6 

T-ReCS 

6530 

1661 1 ® 

2.3 ±0.4 x 10 36 

1.2 -4.6 x 10~ 6 


1 The dust mass for day 6070 has been recalculated from Paper I with the value of « adapted 
to the silicate case. 


2 Time elapsed since the explosion 
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Fig. 1. — (a) SN 1987A seen with T-ReCS at day 6526 in the Si5 narrow band filter (11.7 fxm) 
and (b) at day 6552 in the Qa filter (18.3 /xm); (c) temperature map assuming pure black body 
emission and ? extinction law; (d) Opacity map resulting from the same algorithm. These 
images are smoothed 2 pixels (0.18 arcsec). 
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Fig. 2. — Comparison of the N filter (10 p,m) emission at day 6067 with the Si5 filter (11.7 /vm) 
emission at day 6526: (a) flux calibrated difference of the two images; (b) that difference divided 
by the 2003 image; contours are drawn at [-50,150,250,350,700] percent changes; (c) image with 
contours from the 2003 image; (d) image with contours from the 2005 image. All the images are 
smoothed 2 pixels (0.18 arcsec). 
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Fig. 3. — Overlays of the contour images obtained with ACIS on Jan. 9-13, 2005 at the Chandra 
X-ray Observatory superimposed on the T-ReCS Si5 (a) and Qa (b) image. The T-ReCS image 
has been smoothed 2 pixels (0.18 arcsec). 
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Fig. 4. — Upper: contours from the ATCA image obtained on Jul. 31, 2003 in the 12 mm band 
(16-26 GHz) superimposed on the T-ReCS Si5 (a) and Qa (b) images; lower: ATCA image at 
same frequencies obtained on May 5, 2004 (ATNF Web page) with contours from T-ReCS at 
11.7 pm (c) and 18.3 pm (d). The T-ReCS image have been smoothed 2 pixels (0.18 arcsec). 
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Fig. 10. — The Spitzer IRS spectrum fitted with a silicate dust model. Green circles: Spiteer 
flux measurements. This model has been scaled to fit the T-ReCS flux measurements (red 
squares). Note that the mid-IR emission has been brightening significantly between both sets 
of observations. The Spitzer IRS data unambiguously show that silicate are the major dust 
component. The derived parameters of the fit are: T = I 8 OI 20 K, Mdust = 0.7 — 1.7 x 1 O~ 6 M 0 , 
and Li R = 1.6 ± 0.3 x 10 36 erg s -1 . 







Fig. 12. — Temperature (a) and optical depth r(11.7 fxm) (b) maps derived for the T-ReCS 
observations for a silicate dust composition. 
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Fig. 14. — (a) Montage of the mid-IR 11.7 jim. image of the ER (red) and the HST F625W image 
obtained on day 6502 (yellow); (b) same as (a) with the HST image (yellow) convolved to match 
the resolution of the T-ReCS 11.7 /J,m image (red); (c) montage of the Qa image from T-ReCS 
(red) with the HST F625W image (yellow); (d) same as .(c) with the HST image convolved to 
match the T-ReCS Qa angular resolution. 






Fig. 15. — (a): Image of the ER obtained through the 11.7 gm filter with T-ReCS with contours 
of the HST UVO image obtained on day 6502 superimposed; (b): same with the 11.7 gm image 
deconvolved with the “Multiscale Maximum Entropy Method” from ?. 
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Fig. 16.— Light curves of SN 1987A for the indicated frequencies; the X-ray data points 
before day ~ 4000 are from the ROSAT data (?); ATCA and ROSAT/ACIS fluxes are 
plotted in arbitrary logaritmic scales; ATCA data points are from the ATNF web page 
(http://www.atnf.csiro.au/research/SN1987A/) and ACIS data are from ?. The upper point 
in the Spitzer / IRS data is integrated in the 18.3 pm bandpass, and the lower point is integrated 
in the 11.7 pm bandpass. The point called ‘Ejecta’ derives from the point source near the center 
of the image reported in Paper I. After day 4000 all data points refer to the inner ER. 





iilliiiii 


’" v .' :;;K 




gH 




Fig. 17. — Mid-IR images of SN 1987A from Spitzer, compared to a near-IR image from the 
Hubble Space Telescope. Left to right, top to bottom: (a) Hubble ACS/HRC F814W (range = 
[30, 3.2e4] electrons, pixel = 0.05”); (b) Spitzer IRAC 3.6 fj,m (range=[-0.2, 4] MJy/sr, pixel = 
1.2”); (c) Spitzer IRAC 4.5 /mi (range=[0.2, 5] MJy/sr, pixel = 1.2”); (d) Spitzer IRAC 5.8 /im 
(range=[1.8, 8] MJy/sr, pixel = 1.2”); (e) Spitzer IRAC 8 /mi (range=[2, 20] MJy/sr, pixel = 
1.2”); (f) Spitzer MIPS 24 /mi (range=[18, 45] MJy/sr, pixel = 2.45”). The bright stars in the 
Hubble image NW and SE of SN 1987A (Star 2 and Star 3) appear to be brighter than the SN 
emission at 3.6 /mi, but rapidly fade as the SN brightens at longer wavelengths. The large ring in 
the MIPS image is the first diffraction ring of the PSF, which is visible for any bright unresolved 
source. All pictures are oriented with N upward and E to the left. All span a 35.2” x 34.6” 
field. The Hubble image is a preview image from the Multimission Archive at Space Telecope 
(MAST). The Spitzer images are post-BCD images from the science data archive at the Spitzer 
Science Center. 
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Fig. 18. — IRS slit positions overlaid on the Spitzer IRAC 8 /mi image (range == [0,15] MJy/sr). 
This image spans a 82.8' square field. The slit positions indicated by the long narrow red/yellow 
rectangles are the IRS short-lo (SL) lst/2nd order (7.4-14.2 / 5. 3-8.5 /on) pointings. The short 
narrow blue rectangles indicate the short-hi (SH) pointings. The short wide green rectangles 
indicate the long-hi (LH) pointings. For each IRS module the target is positioned at locations 
1/3 and 2/3 of the distance along the slit length. 






